CLNS 98/1587 
CLEO 98-16 



O 
O 

o 

(N 

00 



> 



o 



Search for Direct CP Violation in S Hyperon Decay 

CLEO Collaboration 

(February 7, 2008) 

Abstract 



Using data collected with the CLEO II detector we have performed a search 
for direct CP violation in the H hyperon system. CP violation gives rise to an 
asymmetry, A, between the parity-violating angular distributions of the decay 
chains H~ — > Avr", A -^ pvr" and H"*" — > Avr"*", A — > pTT~^. In the Standard 



r — . Model, A^ W ^ to 10 ^. If CP violation were found at a higher level it could 

indicate the presence of new physics. We find no evidence for CP violation: 
0\] A= -0.057 ± 0.064 ± 0.039. We also obtain A^ = -0.070 ± 0.064 ± 0.045 and 

§: a= = 0.49l0:°i±0.10. 
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To date CP violation has only been observed in the neutral kaon system [Q and its 
origin remains unknown. It is believed that there is insufficient CP violation in the minimal 
standard model (MSM) to generate the matter-antimatter asymmetry of the universe [0. 
Searches for additional CP violation beyond the MSM may help reconcile this problem. 

In the MSM CP violation effects are due to a single complex phase in the CKM quark 
mixing matrix 0. In the standard phase convention the two matrix elements with large 
phases are V^b and Vm- Because these elements have small magnitudes and involve the third 
generation, CP violation in kaon decays is small. In the kaon system, e is a measure of CP 
violation due to mixing [^ and e' is a measure of CP violation in the decay amplitude, called 
direct CP violation. Theoretical errors are too large to interptret the measured value of e'/e 
as evidence for the MSM mechanism of CP violation or for physics beyond the MSM. In the 
MSM, and its extensions, large CP-violating effects are anticipated in B meson decays. CP 
violation may also be observed by comparing AS = 1 hyperon and antihyperon non-leptonic 
decays |]§J§]. The latter is the subject of this Letter. 

In hyperon and kaon decay gluonic penguin transitions [^ are thought to give rise to the 
relative weak phase difference between particle and antiparticle amplitudes. In kaon decays 
the interference is between the two isospin amplitudes, whereas in S decay it is between the 
5" and P wave amplitudes. Consequently, in S decay the CP violation observable is not 
a rate difference but a difference. A, in the degree of parity violation in charge conjugate 
states. MSM predictions of ^ depend on the values of e, the top quark mass and the hadronic 
matrix element. The dominant uncertainty is due to incomplete knowledge of the hadronic 
matrix element. The Superweak model and models with a very heavy neutral Higgs where 
there are no \AS\ = 1 CP-odd effects predict no CP asymmetries. Other models in which 
lA^"! = 1 CP nonconservation is dominant, predict asymmetries which are on the order of, 
or larger than those in the MSM. Theoretical predictions of ^4= and A\ range from 10~^ to 
10-5 fBHTTI. 



Parity violation occurs in hyperon ^ ^ i decays due to the existence of two orbital 
angular momentum amplitudes of opposite parity. The parity violation observable is an 
asymmetry in the angular decay distribution due to interference between the two amplitudes. 
In the decay S" -^ Ayr" followed by A — *^ P^^', the A is produced with a polarization equal 

to 

(aH + A-P3)A-/5s(AxP3)-7HAx (AxPs) ^^^ 

■T A = ' (1) 

(l + a^A-Ps) 

where Pw is the S polarization, a^, /?= and 7= are the S asymmetry parameters, which 
measure the degree of parity violation in S decay, and A is a unit vector along the A 



momentum in the H rest frame |12]. If the S polarization is unobserved, or if the S is not 
polarized. Equation |I| reduces to Pa = as-A. The angular distribution of the proton from 
the decay of the A is therefore 

dN 

oc 1 + aH«A cos 6'a (2) 



dcos9 



where 6\ is the angle between the proton momentum vector in the A rest frame and A, and 
a\ is the A decay asymmetry which measures the degree of parity violation in A decay. If 



CP is conserved in the decays A ^ pvr and S -^ Att, a^ = — c^a and a= = — ag, respectively. 
Therefore the CP- violating asymmetry parameters, A\ and ^4=, are defined as [^ 

A^ = ^^±^ and A. = ^^±^. (3) 

If CP is conserved in the decay sequence S -^ Avr, A -^ pn, a^aA = ttgc^A- Therefore we 
measure the CP- violating asymmetry parameter WM: 



A = -^—^ ^^ = As + Aa + 0(AaHA«A) -A^ + Aa (4) 

«H«A + «g«A 

where Aa^ a = ^h.a + «h,a- 

The data sample in this study was collected with the CLEO II detector |[1^ at the Cornell 
Electron Storage Ring (CESR). The integrated luminosity consists of 4.83 fb~^ taken at and 
just below the T(45') resonance, corresponding to approximately 5 million e~^e~ — *■ cc and 
1.25 million e+e~ -^ ss events. The latter process is the dominant source of S's in this 
analysis. 

The A is reconstructed in the pn~ decay modeQ We require two oppositely charged 
tracks to originate from a common vertex. The positive track is required to be consistent 
with a proton hypothesis^. The momentum of the A candidate is calculated by extrapolating 
the charged track momenta to the secondary vertex. The invariant mass of A candidates is 
required to be within three standard deviations (3cr = 6.0 MeV/c^) of the known A mass. 
Track combinations which satisfy interpretation as K^ -^ tt+tt^ are rejected. Combinatoric 
background is reduced by requiring the momentum of the A candidates to be greater than 
800 MeV/c. Due to this requirement B decays are not a signifcant source of S's. 

The H~ is reconstructed in the Avr" decay mode. S~ candidates are formed by combining 
each A candidate with a negatively charged track consistent with a pion hypothesis. The S~ 
candidate vertex is formed from the intersection of the A momentum vector and the nega- 
tively charged track. To obtain the S~ momentum, and Att invariant mass, the momentum 
of the charged track is recalculated at the new vertex. The invariant Avr mass is shown in 
Fig. m for all S~ candidates satisfying the selection criteria. 

The Avr mass is fit is to a double Gaussian with widths and ratio of areas fixed from a 
GEANT |T^ based Monte Carlo (MC) simulation of the detector and a first order Chebyshev 



polynomial to describe the combinatorial background. The S — >■ Avr events are generated 
using LUND/JETSET 7.3 [0, in the process e~^e~ —>■ ss —^ EX. The MC momentum 
distributions and production angles of all particles in the decay chain are in good agreement 
with observation. We find 8434 ± 109 events consistent with H -^ Avr. The mean of the Att 
invariant mass distribution is in agreement with the MC simulation. 



^Charge conjugation is implied in this paper except were explicitly noted 

^Hadronic particles are identified by requiring specific ionization energy loss measurements 
(dE/dx), combined with time-of- flight (TOF) information when available. The two measurements 
are combined into a joint probability for the particle to be a pion, a kaon or a proton. A charged 
track is defined to be consistent with a particle hypothesis if its probability is greater than 0.003. 
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FIG. 1. Att invariant mass distribution. 



The cos6'a resolution function has both Gaussian and symmetric non-Gaussian compo- 
nents. The resolution in cos6'a, aav, defined to be the average of the rms variance of the two 
components, weighted by their relative normalizations is aav = 0.030. 

The product of the decay asymmetry parameters, a^a^, is measured using a one dimen- 
sional unbinned maximum likelihood fit to the cos 6'a distribution given in Equation |], in a 
manner similar to [|18|. This technique enables a likelihood fit to be performed to variables 
modified by experimental acceptance and resolution. The probability function of the signal 
Ts is determined by generating one high statistics MC sample of S~ -^ Avr" at fixed asttA- 
By suitable weighting of the accepted MC events a likelihood is evaluated for each data event 
for trial values of a^a^, and a fit performed. The probability for each event is determined 
by sampling Ts using an interval centered on each data point. 

Background is incorporated directly into the fit by constructing the log-likelihood func- 
tion: ln£ = J2iLi ^^(PsTs + Pb^b) where N is the number of events in the signal region, 
defined to be within 7.5 Mev/c^ (iSa) of the S mass, and Ps and Pb are the probabilities 
that events are signal and background respectively. The probability distribution of back- 
ground in the signal region, F^, is determined from 15 MeV/c^ wide (Avr) mass sidebands 
above and below the signal region. 

To ensure that there is no gross bias in the analysis technique we assume CP conservation 
and determine the decay asymmetry in S~ — >■ Avr", combining S"*" and S~ events. We find 
a= gctA^A = —0.291 ± 0.021 where the error is statistical. The result is in good agreement 
with the world average asttA = —0.293 ± 0.007 [0 . 

We have considered the following sources of systematic uncertainties and give our esti- 
mate of their magnitude in parentheses. The error associated with finite MC statistics is 
estimated by varying the size of the MC sample used in the fit (1.6%). The uncertainty 
from varying the A selection criteria and particle identification is (—6.3%, +2.1%). The un- 
certainty associated with MC modeling of slow pions from S and A decay is obtained by 
varying the slow pion reconstruction efficiency according to our understanding of the CLEO 
II detector (—0.0%, +0.7%). The error due to finite interval size is determined by varying 
the size of the interval for both signal and background (—0.5%, +4.1%). The uncertainty 
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FIG. 2. The sideband subtracted angular distribution of the proton in the A rest frame with 
respect to the direction of the A in the H~ rest frame in H~ -^ Avr^, A — > pn^ for data (points with 
error bars) and projection of the fit (Une). Both charge conjugate states are included. 



arising from incomplete knowledge of the background shape in the sidebands is evaluated by 
varying the size of the mass sideband region (—1.2%, +0.9%). A binned fit is also performed 
and the result is consistent with the maximum likelihood fit. This measurement is insensitive 
to production polarization, Pw, and no systematic error has been included from this source 
(see below). Combining all sources of systematic errors in quadrature the total systematic 
error is (-6.6%, +5.0%) and aH,gaA,A = -0.291 ± 0.021(stat)+|!;;}5(syst). 

To search for CP violation the S candidates are sorted into charge conjugate states. There 
are 4204 + 75 S~ decays and 4200 + 75 S+ decays. Using the same fitting procedure as above, 
and one MC sample for each charge conjugate state, we measure a^aA = — 0.275 + 0. OSOIqIoJs 
and a^a^ = —0.308 + 0.030lo;oi6 respectively. Data and fit projections are shown in Fig. |[ 
From these results we calculate A = 0.057 + 0.064 where the statistical error is estimated 
using the method in 



14 



Most systematic errors associated with S~ reconstruction apply 
equally to S+ and cancel when determining the asymmetry A. If there is any S production 
polarization, it is required to be normal to the production plane. The angular distribution 
with respect to this plane will have opposite slopes for S~ and H'^. However both the 
production plane and detector acceptance are uniformly distributed in the azimuthal angle 
0. Therefore this result is insensitive to production polarization. If there is any difference 
between the relative angular efficiencies of S~ and S+ that is not modeled accurately by MC, 
this would lead to a biased value of A. A momentum dependent asymmetry in detection 
efficiencies for negative and positive soft pions could, in principle, exist due to differing 
hadronic interaction cross sections of vr" and vr^ in the 3.5 cm radius CLEO II beryllium 
beam pipe. We assume this effect to be absent for two reasons. (1) In this analysis, typical 
decay lengths for S and A are 6 cm and 11 cm respectively, therefore most tracks originate 
outside the beam pipe. (2) There is no evidence for a momentum dependent asymmetry in 



the CLEO II detection efficiency for tracks that originate inside the beampipe |T9 |. 

The A7r~(A7r+) mass sidebands consist mostly of real A's paired with random pions. 
The efficiency corrected angular distribution of the sidebands in each charge conjugate state 
will be the same if the MC accurately models the relative angular efficiency of slow vr"*" 
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FIG. 3. The sideband subtracted angular distribution of the p in the A rest frame with respect 
to the direction of the A in the H~ rest frame in H^ — > Avr^, A — > pn^ for data (triangles with error 
bars) and projection of the fit (dashed line). The sideband subtracted angular distribution of the 
p in H"*" — > Avr"'", A — > pvr" in data (squares with error bars) and projection of the fit (solid line). 

and 7i~ . To study this precisely, a high statistics sideband sample is obtained by relaxing 
the selection criteria and increasing the width of the sideband regions. The fit, performed 
on 12365 ± 111(12353 ± 111) events in the A7r^(A7r+) sideband regions, yields a=aA = 
0.062 ± 0.016 in the S" sideband and agOA = 0.077 ± 0.016 in the H+ sideband. The 
difference between these two values —0.015 ± 0.023 indicates that there is no statistically 
significant discrepancy in the modeling of the H~ and S"*" detector angular acceptance in MC 
compared to data. The error on the difference is taken as the systematic error on the CP 
asymmetry parameter measurement. The result is 



A 



-0.057 ±0.064 ±0.039 



(5) 



There are three measurements of A/^ |2y] [^] |£2|. Using the most precise: A/, = ±0.013 ± 
0.022 p^ we obtain A^ = —0.070 ±0.064 ±0.045. Where the second error is systematic and 



includes the error on A^. 

It is possible to obtain a more precise value for A by combining our value for aga^ 
with the PDG value of a^aA to obtain A = -0.025lg:i]^^ and A^ = -0.038t[]:J]^^ where the 
statistical and systematic errors have been combined. 

Finally, as this measurement of aga^ is significantly more precise than the only published 



value ^], by combining with: a^ = -0.63±0.13 0, we obtain ag = 0.49to:o5±0.10, where 



the first error is the combined statistical and systematic error from this study and the second 
error is from a^. 

In conclusion, from a sample of approximately 4,000 S decays and similar number of S, 
we have measured agOA = -0.308 ± OmOt^fii and obtained ag = 0.49l[J:[]^ ± 0.10. We 
have searched for direct CP violation by measuring the asymmetry parameter A. We find no 



evidence for CP violation: A = -0.057 ± 0.064 ± 0.039 and A^ = -0.070 ± 0.064 ± 0.0450. 
We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss National 
Science Foundation, and the Alexander von Humboldt Stiftung. 
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As this paper was about to be submitted, we became aware of K.B. Luk et al. hep-ex/0007030| , 



July 2000. This paper measures aga^, A and ^4= with greater precision than the values reported 
here. The two analyses are consistent but employ different experimental techniques. 
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